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Abstract. We present the analysis of FUSE observations of two Herbig AeBe stars known to harbor young 
circumstellar disks: HD 100546 and HD 163296. In both cases we detect absorption lines from warm and dense 
H2. The thermalization of the rotational levels up to J ~ 4 allows evaluation of the temperature, density and 
typical size of the absorbing layer. These quantities are consistent with absorption of the light of the central star 
by a thin layer of a circumstellar disk seen at an intermediate inclination. 
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1. Introduction 

HD 100546 and HD 163296 are two Herbig AeBe stars 
harboring young circumstellar disks. Herbig AeBe stars 
are intermediate-mass (2-10 Mq) young emission-line 
stars still in the pre-main sequence phase (Herbig [i960). 
Material surrounding Herbig stars is thus considered as 
the constituent of protoplanetary disks, similar to the 
disks seen around the less massive T-Tauri stars. These 
disks are the material from which (and the places where) 
planets are supposed to form according to the standard 
model of planetary formation (Lissauer|1993). 

HD 100546 (B9Vne, d=103t£pc, t>10Myr; van den 
Ancker et al. H55S)) is older than HD 163296 (AlVe, 
d=122±^pc, i=4±| 5 Myr; van den Ancker et al. 11555)1 . 
However, both stars present very similar characteristics. 
In the two cases, ISO detected silicate features resem- 
bling cometary materials (Malfait et al. 1998. van den 
Ancker et al. 125551 Bouwman et al. 125551 TWO} . Disks 
have recently been imaged around both stars (Pantin et 
al. 125551 Grady et al. 125551 1255T1 Augereau et al. !255T)l . 
The images show dusty disks at intermediate inclination 
(~60°). Absorption spectroscopy also revealed the pres- 
ence of circumstellar gas around HD 100546, including 
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spectral variabilities which are interpreted as accreting 
material (Grady et al. H555I IT557)l . 

In the case of HD 163296, a CO disk of about 310 AU 
in semi-major axis has been detected (Mannings & 
Sargent 1997J). H2 infrared emission lines have also been 
reported by Thi et al. ( 2001.) at 17 and 28/im. From their 
ISO observations and assuming thermal equilibrium up to 
J = 3, they derive a total H 2 mass of 0.4(±0.2) x 10 -3 M , 
and the ratio of the two emission lines (v = 0, J=2^0 
and J=3— >1) gives a temperature of 220(±22) K. However, 
this detection has been challenged by more recent ground 
based observation at 17/im which shows no detection with 
three times better sensitivity (Richter et al. P002 ) . 

These controversial observations raise the question of 
the molecular content of such disks, where most of the 
mass is in the form of molecules, and particularly H2. 
The formation of giant planets requires a large reservoir 
of molecular gas. Moreover some of these planets migrate 
close to their parents where they are observed. This mi- 
gration also needs a massive disk to allow the angular 
momentum exchange between the migrating planet and 
the disk. The H2 content is thus a key ingredient in the 
recipe for the formation of giant planets. But H2 is a sy- 
mctrical molecule, and infrared emission by quadrupole 
rotational transitions is very inefficient. In contrast, when 
seen against a UV bright source, far-UV absorption lines 
due to electronic transitions allow sensitive observation 
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of H 2 even at low temperatures. FUSE offers a unique 
opportunity to scrutinize in detail the H 2 content of pro- 
toplanetary disks around Herbig stars. 

2. Observation and data analysis 

HD 100546 was observed with FUSE through the LWRS 
aperture (30"x30") for a total time of 5.8 hours on 
March 26, 2000 (Program P1190303) and March 3, 2002 
(Program P2190401) (for an overview of FUSE, see Moos 
et al. I2UUUI and Sahnow et al. I2UUU|| . HD 163296 was ob- 
served twice, on April 27, 2001 and April 29, 2001 for a to- 
tal time of 8.9 hours (Programs P2190601 and Q2190101, 
respectively). The data of both targets were reprocessed 
with the version 2.0.5 of the CALFUSE pipeline. The out- 
put of the pipeline is a total of 5 and 10 sub-exposures 
for HD 100546 and HD 163296, respectively. The sub- 
exposures have been aligned and coadded resulting in a 
set of four independent spectra, one for each FUSE chan- 
nel (2 LiF spectra and 2 SiC spectra). 

The version 2.0.5 of the CALFUSE pipeline is known to 
slightly over-estimate the tabulated errors on each pixel. 
We compared the data used in the present work with 
data of the same observations but obtained with the ver- 
sion 2.2.1 of the CALFUSE pipeline; we conclude that the 
error bars given below are not significantly affected by 
the improvement in the error propagation of the different 
pipelines. 

Apart from the observation of the molecular hydro- 
gen, the spectra are rich in emission and absorption lines 
from atomic and ionic species. For instance, a large num- 
ber of Fe 11 lines from the ground level as well as from 
excited levels are clearly detected. As already observed 
in AB Aur and (3 Pic (Roberge et al. 120011 Deleuil et 
al. 1200 ljl . bright emission from Cm and Ovi are also de- 
tected in the two spectra. All these features indicate the 
presence of circumstellar material, chromospheric activity 
and/or accretion (Bouret et al. I2002f) . A detailed analysis 
of the whole FUSE spectrum will be made in a forth- 
coming paper (Deleuil et al. 2003). Here we focus on the 
analysis of the H 2 lines, probing the molecular portion of 
the disk. 

In the far-UV, HD 100546 is brighter than HD 163296 
and the observed spectrum has a better S /N ratio allowing 
detection of H 2 in its pure rotational level (v = 0) up to 
J = 9 (a line of v = 0, J = 10 is also marginally detected 
at 1058.6 A) and in its first vibrational level (v — 1) up 
to J = 5 (Fig. P). In the HD 163296 spectrum, H 2 lines 
are detected up to v = 0, J = 4. Higher J-levels are 
beyond the detection limit. The data analysis and pro- 
file fitting has been done using the Owens code kindly 
made available to us by Dr. M. Lemoine (see for exam- 
ple Lemoine et al. 12303 and Hebrard et al. l2Ho^jl . For the 
electronic transitions, we used the wavelengths and oscilla- 
tor strengths tabulated by Abgrall et al. ( 1993a and 1993b 
for the Lyman and the Werner system, respectively), and 
the inverses of the total radiative lifetimes tabulated by 
Abgrall et al. plOTT)) . 
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Fig. 1. Sample of H 2 lines detected toward HD 100546. 
The histogram shows the data. The final fit to the data is 
shown (thick line) together with the theoretical spectrum 
without convolution with the instrumental line spread 
function (thin line) . This shows that the fitted lines are not 
saturated. In the selected sample, some lines are blended 
together. Consequently, a profile fitting is needed to use 
the information included in the data but which can not 
be translated into an equivalent width. 

With low H 2 column densities in HD 100546 (N < 
10 16 cm -2 for each J-level), we can choose to fit only 
the unsaturated lines having low oscilator strength (/ <, 
10 -2 ). This allows us to avoid systematic errors in the es- 
timates of column densities which could be included by 
the fit of saturated lines with uncertain instrumental line 
spread function. Note, however, that additional fits have 
been performed with saturated lines (/ *b 10~ 2 ). They 
give very similar column densities and allow the determi- 
nation of the intrinsic line width (b). In HD 163296 the 
column densities are larger, and all the H 2 lines are sat- 
urated. Therefore the error bars on the column densities 
are also larger than in HD 100546. An example of the fit of 
two H 2 lines ( J — 3 and J — 4) which are superimposed 
on the O VI stellar emission line is given in Fig. [21 

The estimated column densities are tabulated with 2a 
error bars (Table We obtain the total column densities 



and 



N H2 (HD 100546) 



A^h 2 (HD 163296) 



2.8 x 10 16 cm" 2 



1.4 x 10 18 cm" 2 . 



The error bars are estimated by the classical method 
of the A% 2 increase of the x 2 of the fit (see Hebrard et 
al. 1)2002(1 for a full discussion on the fit method and er- 
ror estimates with owens profile fitting). These error bars 
include the uncertainties in the continuum fits, intrinsic 
line widths (b) and the instrumental line spread func- 
tion. We note that the observation of different lines with 
different oscillator strengths for the same species allows 
us to constrain all these quantities. In particular the fit 
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Fig. 2. Two H2 absorption lines detected over one emis- 
sion line of the O vi doublet toward HD 163296. Here the 
theoretical spectrum without convolution with the instru- 
mental line spread function (thin line as in Fig. 1) shows 
that these H 2 absorption lines are saturated. 



Fig. 3. Excitation diagram of the H2 lines detected to- 
ward HD 100546. The column densities of the rotational 
levels with J < 4 are consistent with a single temperature 
of T ~ 740 K. This shows that the H2 volume density 
must be higher than the critical value of nn 2 ~ 10 4 cm -3 . 



to saturated lines constrains the line widths and the in- 
strumental line spread function. This last one is allowed 
to slowly vary with the wavelength. The results are ob- 
tained from a final self-consistent fit to all the data in 



which the continuum, the instrumental line spread func- 
tion, and physical parameters of the absorption lines are 
free parameters and estimated by the determination of 
the best x 2 m t ne parameters' space. For the intrinsic 
line widths, we obtain &hdioo546 = 3.4+Q?km s _1 and 



Table 1. H 2 column density toward HD 100546 and &hd 163296 — 2. 2_ ; 7 km s 

3. Results 

3.1. Excitation diagrams 
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+0.4 
1.0 
+0.2 
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13.87 
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The error bars are 2— a error bars. 



The excitation diagrams corresponding to the derived 
H 2 column densities are presented in Fig. [3] and 01 for 
HD 100546 and HD 163296, respectively. For the v = 1 
levels observed toward HD 100546, we see that the column 
densities are above a simple extrapolation of the v = lev- 
els as plotted on the excitation plot (Fig. 0J. This effect 
is also observed in the vibrationally excited interstellar 
H 2 detected toward HD 37903 (Meyer et al. l2DDT)l . where 
the high energy levels are not thermally populated. Note 
also that the level degeneracies (g) of the v = 1 levels 
are smaller than the degeneracies of the v = levels with 
similar energy, hence in the same region of the plot. The 
plotted values of N/g are thus larger for these v = 1 levels 
than for the nearby v = levels (Fig. 

3.2. Turbulence 

The simultaneous determination of temperature and line 
broadening can be used to constrain the turbulence of the 
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Fig. 4. Excitation diagram of the H2 lines detected to- 
ward HD 163296. The column densities of all the de- 
tected J-levels are consistent with a single temperature of 
T ~ 430 K. This shows that the H2 volume density must 
be higher than the critical value of h,h 2 ~ 8 x 10 4 cm~ 3 . 

gas. The observed line width b is a combination of the 
thermal broadening, the intrinsic turbulence Wturb , and the 
projected radial component of the gas motion seen along 
the star diameter. With the knowledge of the temperature, 
and considering that the motion can only increase the line 
width, we can estimate a lower limit on the turbulence 
velocity of the gas. We have 

[2kT ~ 

— A/ turb 
V MH 2 

where /in 2 is mass of molecular hydrogen. This can be 
written as 

f^y<(^y- 8-3ixio-^. 

Using the temperature given above and the 2a up- 
per limits on the measured line widths (&jjd 100546 — 
4.2kms~ 1 and &HD163296 — 2.7kms _1 ), we find that 

Wturb,HD100546<3.4kmS _1 and U t urb,HD 163296<l-9kmS _1 . 

3.3. Nature and origin of the H2 gas 

Many characteristics show that the detected H 2 is not in- 
terstellar but circumstcllar; this favours the interpretation 
that the detected lines arise from protoplanetary material. 
First, in both cases the radial velocity of H2 is similar to 
the radial velocity of the atomic lines due to the circum- 
stellar gas. Lines from exited levels of Fen (like Fen* and 



Fell**, Deleuil et al.[2003) cannot be due to the interstel- 
lar medium. They must be linked to dense circumstellar 
material. The similarity of their radial velocity with the 
radial velocity of H2 is a first clue that the H2 is linked 
with the circumstellar matter. 

Importantly, the excitation diagrams show that the H2 
is thermalized up to about J — 4, with high temperatures. 
Hot H2 has also been observed around the T-Tauri star 
TW Hya (T=3000 K, Herczeg et al. EMty . Here we obtain 

T w 740 ± 30 K for HD 100546 

and T « 430 ± 20 K for HD 163296. 

Seen in absorption, it is very unusual to observe such 
a thermalization at high temperature up to J = 4 (to 
our knowledge, this is the first case). This shows that 
the absorbing medium has very particular physical condi- 
tions, different from what is commonly seen in the diffuse 
interstellar medium. Indeed this requires a line of sight 
which presents simultaneously high density for the ther- 
malization and low thickness to have column densities low 
enough for a detection of the far-UV stellar continuum. 
In the diffuse and transluscent interstellar media, levels 
higher than J — 1 or J — 2 are populated by mechanisms 
other than the collisional pumping. In these low density 
media, the collisional pumping is less efficient than the 
UV-pumping and mecanical processes such as C shoes or 
turbulence. Here levels up to J = 4 show a dominant col- 
lisional pumping. This is a clue of a high density medium. 

Note here that for the levels higher than J = 4 as 
detected in HD 100546, it is not possible to discriminate 
between a collisional pumping within a higher tempera- 
ture component of the absorber and other mecanical or 
radiative processes. 

The observed high temperatures give additional clues 
that the detected H2 is circumstellar and close to the ex- 
citing stars. Using the numerical model developed by Le 
Bourlot et al. (1993), we find that the incident radiation 
on the H2 must be ~ 10 4 times larger than the mean galac- 
tic UV radiation field. This shows that the observed H2 
is circumstellar and really close to its central star. These 
temperatures correspond to ~1.5 AU from HD 100546 and 
-4 AU from HD 163296. 

Most importantly, the thermalization allows us to es- 
timate the H2 volume density (nn 2 ) . Using the critical 
density given by Le Bourlot et al. I|1999|) . we find : 

n H2 (HD 100546) > 10 4 cnU 3 

and n U2 (HD 163296) £ 8 x 10 4 cm -3 . 

Assuming these volume densities, we can estimate 
the typical size of the absorbing layer responsible for 
the observed absorptions. Note that with low column 
densities and large volume densities, the absorbing cloud 
must be thin. The ratios of the H2 column densities to 
the volume densities give a typical thickness of 550.2 AU 
for HD 100546 and £1.2 AU for HD 163296. 
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All the observed and estimated quantities are consis- 
tent with absorption of the star light by a thin outer layer 
of a disk seen at an intermediate inclination. Because the 
disks of HD 100546 and HD 163296 are both inclined by 
about 60° , it is likely that we observed a line of sight graz- 
ing the disk surface. In any case, this shows that the H2 
gas is not confined to a very flat disk. Alternatively, there 
is still a possibility that there might be shocked molecu- 
lar gas in the vicinity of these stars, due to the combi- 
nation of a bipolar jet (as observed around HD 163296) 
and a circumstellar envelope. However it is not possible 
to discriminate between these two possibilities with solely 
spectroscopic observation with no spatial information. 

4. Conclusions 

We observed the molecular hydrogen toward two Herbig 
stars surrounded by young protoplanetary disks. The ob- 
servation of an extremely large number of rotational and 
vibrational levels allows for the first time the determina- 
tion of the H2 physical conditions within circumstellar gas 
around two HAeBe stars. With this detailed view of the 
population of the H 2 excited levels, we now need a full 
modeling in which all the cooling and heating terms must 
be taken into account. Such a modeling is under devel- 
opment. This opens a new window into the physical con- 
ditions operating in the protoplanetary disks surrounding 
massive stars. 

The present observations of material surrounding stars 
with typical age between 10 6 and 10 7 years give new 
constrains on the relation between the radiation and the 
molecular protoplanetary disk in which planets are sup- 
posed to form and migrate (Terquem et al. I2()0()|l . This 
gives information on the interface between the dense part 
of the disks and the interstellar medium, and can help to 
better understand the late stages of the planetary forma- 
tion. The radiation is indeed a key element in the evapo- 
ration of the protoplanetary disks whose lifetime seems to 
be limited (Zuckerman et al. 1995) although the process 
responsible for the clearing of the disks remains unclear. 
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